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^ Abstract 

in 

m . 

psj , The sinh-Gordon model with integrable boundary conditions is considered in low order 

O ■ perturbation theory. It is pointed out that results obtained by Ghoshal for the sine- 

' 

\ Gordon breather reflection factors suggest an interesting dual relationship between models 
with different boundary conditions. Ghoshal's formula for the lightest breather is checked 
perturbatively to 0((3 2 ) in the special set of cases in which the </> — > —<j) symmetry is 
maintained. It is noted that the parametrisation of the boundary potential which is natural 
for the semi-classical approximation also provides a good parametrisation at the 'free- 
fermion' point. 

W. 

1. Introduction 

In recent years, some progress has been made concerning the question of integrability 
for two-dimensional field theories restricted to a half-line [fT|-|10||, or to an interval. (For a 
recent review of some aspects of this, see In this article it is intended to concentrate 

on the simplest example, other than free field theory, and to estimate using perturbation 
theory the reflection factors corresponding to those boundary conditions which preserve 
integrability. Thus, the model to be studied is the sinh-Gordon model restricted to the 
region x < 0. There are several reasons for doing this. Firstly, Ghoshal || has suggested 
a formula for the reflection factors for the breather states of the sine- Gordon model using 
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the reflection bootstrap equations and the soliton reflection factors suggested earlier by 
Ghoshal and Zamolodchikov It is plausible to suppose that the reflection factor for 
the lightest breather state, analytically continued to imaginary coupling, should provide 
a sensible hypothesis for the reflection factors for the sinh-Gordon particle. On the other 
hand, the reflection factors of the sinh-Gordon model are accessible perturbatively and so 
a start can be made towards checking Ghoshal's formula. Unfortunately, the relationship 
between the two classical parameters which may be introduced via the boundary poten- 
tial and the quantum reflection factors is unclear except in certain cases. Therefore, the 
perturbative calculation is potentially a useful way to suggest the missing connection. In 
those cases where the relationship between classical and quantum parameters is already 
known, perturbation theory should provide a verification. The second reason concerns 
weak-strong coupling duality. 

The S-matrix describing the elastic scattering of a pair of sinh-Gordon particles is known 
||12j| to have the form 

S{e) --jmhB) (L1) 

where 

1 3 2 

B = - ^r— (1.2) 

2tt 1 + /? 2 /4tt V ! 

and 

(x) = — 2 + 4 . (1.3) 
V ' sinh(f -*f) V ; 

The rapidity difference of the two particles is represented by 0. This is clearly invariant 
under the transformation 

Alt 

^j, d.4) 

and this invariance is referred to as the weak-strong coupling duality. This form of duality 
is shared by all the other affine Toda models based on the ade series of data [jl2| , |13|| , and 
a generalised form of it is enjoyed by those models based on the data derived from the 
other (ie not simply-laced) affine diagrams []14]JI5[] . In some earlier works concerning the 
reflection factors |3|,|| it was suggested that they should enjoy the same duality property. 
However, this is unlikely to be the the following argument demonstrates. 

When there is a boundary condition at x = 0, provided it is chosen to maintain integra- 
bility, then it is expected that a single particle approaching the boundary will be elastically 
reflected from it so that 

|0>out=£(0)|-0>in, 
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where 9 is the rapidity of the particle. In the sinh-Gordon model there is only one such 
particle. Ghoshal's formula for the sinh-Gordon reflection factors, after translation to 
the above notation, has the form: 

K(ff\ = (l)(2-B/2) (l + B/2) 

1 ] (1 - E{(3)) (1 + E{(3)) (1 - (1 + F{(3)) ' 1 ' ' 

where E and F are two functions of f3 which also depend on the boundary parameters, 
but which are not yet fully determined in terms of those parameters.0 However, for the 
Neumann boundary condition Ghoshal requires 

F = 0, E = l-B/2 (1.6) 

leading to a reflection factor /Cat, given by: 

(1 + B/2) 

Kn ~ (B/2) (1) • (L7) 
Clearly, this is not self-dual; rather its dual partner is K.* N , given by: 

Kn = (1-B/2)(1)" (L8) 

What about a verification of (|1.6|) perturbatively? In effect this has been provided 
implicitly by Kim [16, who has developed perturbation theory for affine Toda models 



satisfying the Neumann boundary condition, and has provided detailed calculations up 
to 0((3 2 ). In other words, the 0(f3 2 ) check can be made merely by expanding (|1.7|) and 
comparing with the formulae given by Kim |16|]. Explicitly, from ( |1 . 7| ) , 

1 1 



if3 2 smh6 
/Cjv ~ 1 



cosh 9 — 1 cosh 9 



(1.9) 



which is in perfect agreement.! 

To discover which boundary condition is dual to the Neumann condition it will be enough 
to consider the classical limit of ( |1.8|) . However, first it is necessary to consider the classical 
limit of (|L5|). That is, taking p -»• 0, 

(l) 2 

= ~ (1 - E(0)) (1 + E(0)) (1 - F(0)) (1 + F(0)) ' (L10) 



1 In Ghoshal's notation E = Bt]/tt, F = iB-d/n. 

9 (1) 

Kim did not comment on this agreement. He chose to assume instead that the a\ reflection 
factor should be self-dual and proposed an alternative formula for /Cat, disagreeing with Ghoshal, 
but nevertheless in agreement with perturbation theory to 0(/3 2 ). 



3 



Fortunately, this may be calculated independently and the result was already reported 
some time ago ||. 

In order to state the result of the 'classical' calculation further details of the model are 
required. Its equation of motion and boundary condition are: 

d ^ = _/y| L^H _ e -V2 P A x < 

r (1-11) 

g = -^(a ie ^-a e-^) x = 0, 

where the a/2 is a conventional normalisation in Toda theory. The two constants o"o and 
o"i are essentially free and represent the degrees of freedom permitted at the boundary 
0. The constraints on the boundary parameters are discussed in || and by Fujii and 
Sasaki [ffSJ. The boundary condition is required to have the given form as a consequence 
of maintaining integrability on the half-line. It is also convenient to write 

Oi = cos diix. (1-12) 

Then, the classical reflection factor is determined by first finding the lowest energy static 
solution to fll.llp (the 'vacuum configuration') and then solving the linearised scattering 
problem in this static background. Further details of this are given below and the result 
for the reflection factor calculated in this approximation (and for |a^| < 1) is: 

(D 2 

Ufl =-7 77 -T, 77 7- (1-13) 

(1 — oo — oi) (1 + a + ax) (1 - a + ai) (1 + a - aij 
The formula (|1.13| ) clearly has the same form as the classical limit of Ghoshal's formula, 
and is unity when ao = a± = 1/2 (the Neumann condition), in agreement with the classical 
limit of )Cn, eqfll.T]). 

On the other hand, the classical limit of KL* N is not unity. Rather, it is — 1/(1) 2 , and 
is in agreement with ( |1.13| ) for the choice ao = a\ = 0. In other words, it is tempting to 
conclude that the two boundary conditions 

|^ = and ^ = -^( e WV2_ e - W V2\ at x = ^ (L14) 
ox ox p V / 

are dual, in the sense that their reflection factors are transformed into each other under 

the transformation ( |1.4j ). 

It should now be clear that it is the functional dependence of E and F on the three 

parameters ao, ai and (3 which must be found together with a proper understanding of 

which pairs of boundary data are dual to each other, and which (if any) are self-dual. The 

rest of the paper is concerned with these questions. 
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2. The propagator in a general background 

To calculate the propagator it will be necessary to consider the linear perturbation 
around the static background solution to ( |1 . 1 1|) . The static solutions have the form 

, r- 1-1- P 2(x-X ) 

1 _ e 2(x-x Q ) v > 

where the parameter xq is determined by the boundary condition and is given by: 



cothx /I±^°. (2.2) 
V l + ox 

Note that provided a > c?\ it is guaranteed that xq > 0; if on the other hand o\ > o~o it 
would be necessary to take the background provided by the inverse of the right hand side 
of ( |2.1[ ). In the calculations which follow it will be assumed o~o > o\ for defmiteness. 

Linearising ( |1.11| ) around the static background leads to a pair of equations for the first 
order correction to d>: 



d 2 ^ +4(1+ . ~ ] 0i = x < 

sinn 2{x — Xq) 



dx V V 1 + °i V 1 + Co 



-+(-n/^ + -oJ^)^i = x = 0. 



(2-3) 



Fortunately, these equations may be solved exactly and the required Green's function 
can be written down explicitly once the eigenfunctions of the second order differential 
operator in the first of eqs( |2.3| ) have been determined. It is convenient to write (pk,cu to 
denote the eigenfunction corresponding to the eigenvalue uj 2 — k 2 — 4, and then 

(pk,u = ie~ iujt r(k) (F(k, x)e ikx + F(-k, x)e~ ikx ) , (2.4) 

where r(k) is a real, even function of k (which will be determined by normalising the 
Green's function properly), and where F(k,x) is given by: 

F(k, x) = P{k) (ik - 2 coth 2(x - x )) , 

with 



P(k) = {ik) 2 - 2zfcVl + <W1 + o-i + 2(a + <ti). 



The classical reflection factor ( |1.13|) follows immediately from (|2.4| ) on setting 



2cosh^, k = 2sinh#, calculating F(—k,—oo)/F(k,—oo), and using the convenient 



parametrisation ( |1.12| ). 



The Green's function, or configuration space propagator, constructed from these eigen- 
functions is not hard to find. It is: 



G(x, t- x',t') = i j ' ^ j " ^ g(u, k; x, t; x\ t') (2.5) 

where 

, (f _ t/) f(k, x)f(-k, x ')e ik ^- x ') + K f{-k, x)f{-k, x ') e - ik ^ +x ') 

(2.6) 



G(u,k;x,t;x',t') =e , , . 

— k z — 4 + «e 



... . ik — 2 coth2(x — xq) 

f(k,x) = : , 

JK ' ' ik + 2 

and 

_ P(-fc) zfc-2 
°~ P(fc) zfc + 2' 

When = 2sinh6> this expression for /Co reduces to the previous one, eq( |1.13|) . To check 
that (|2.5|) is normalised correctly it is enough to note that it is in the limit x,x' — > — oo. 



If o"o = <J\ = cr, eq (|2.5|) simplifies because in that case f(k, x) = 1 and the factor JCq 
reduces to 

__ ifc + 2cr 
ik — 2a 

In other words, in this special set of cases, 

p ik(x — x') _i_ t~ p —ik{x-\-x') 

fc; x, t; = e --(*-*') ! t 2 ( 2 - 7 ) 

In this situation, the static background is simply <po = 0. 



In this article a pragmatic approach will be taken [[16]]. The classical reflection factor 
is an integral part of the free field propagator calculated within the classical background. 
The full two-point function may then be calculated perturbatively (once subtractions are 
made to remove infinities), and provides a working definition of the quantum reflection 
factor in the sense that it will be the coefficient of e~ lk ( x+x ) in the residue of the on-shell 
pole in the asymptotic region x, x' — ► — oo. 
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3. Low order perturbation theory 

Since perturbation theory will be developed around an x-dependent background field 
configuration, it is possible to use a more-or-less standard Feynman expansion in configu- 
ration space, using the above propagator but recognising that the 'vertices' will be position 
dependent. One consequence of this is that although the sinh-Gordon theory on the whole 
line has only vertices at which an even number of lines join, the theory on the half-line must 
also, in general, contain vertices at which an odd number of lines meet. In other words, the 
symmetry (p — > — </>, enjoyed by the bulk theory, is broken by most boundary conditions. 
The exceptions to this are precisely the boundary conditions for which o"o = ay. In those 
cases, the classical background is 4>q = and there will be no additional vertices compared 
with the bulk theory other than those arising specifically from the boundary term itself. 
Moreover, the boundary terms will only generate vertices at which an even number of lines 
meet. In this article, it is intended only to perform detailed calculations in the special case 
o"o = <J\ = a. The more general situation will be treated elsewhere since the computations 
are more lengthy and intricate. 

To 0((3 2 ), there are two contributions which need to be calculated and each of them is 
a single loop diagram. The first comes from the boundary and is described by: 

/oo 
dt"G(x, t- 0, t")G(0, t"; 0, t")G(0, t"; x', if), (3.1) 
- oo 

where G\ is a combinatorial factor. The second comes from the bulk potential and has the 
form 

/oo pO 
dt" / dx" G{x,t;x" ,t")G(x" \t";x" ,t")G{x" ,t";x' ,t'), 
-oo J oo 

(3.2) 

where C2 is another combinatorial factor. The relative factor of 8 between the two ex- 
pressions arises from the different factors of a/2 in the bulk and boundary potential terms. 
The expression ( |3.1| ) is relatively simple to analyse and it will be considered first. 

The integral over t" in expression (]3.1|) gives a delta function allowing the u/, say, integral 
to be done leaving the momentum integrals for all three propagators, and the frequency 
integral for the loop (middle propagator) to be performed. 

The middle propagator is clearly divergent and a minimal subtraction will be made to 
yield a finite part. That is, the logarithmically divergent integral 

[ dw" r dk" i f ik" + 2a\ 

J ~2k J ~2^iv" 2 - k" 2 - 4 + ie \ + ik" -2a) 
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will be replaced by the finite part 

du" f dk" i ( 4a 



2n J 2n u" 2 - k" 2 -4 + ie \ik" - 2a 
which can be evaluated to obtain first 

1 f dk" 1 / 4a 



2 J 2tt y/ kl/ 2 + 4 \ik"-2a 

and then finally, using ( |1.12| ), to obtain 

( a -n cosan a i\ 

— (a mod 1) . (<j.o) 

sina7r 

The expression ( |3.3| ) is invariant under the transformation a — > a + 1 . 

The momenta of the other two propagators may be integrated out by closing the contours 
in the upper half plane (first setting k — > —k in the first term of the first propagator). 
Provided a > 0, the reflection factor term does not contribute a pole. On the other 
hand, if a < 0, there is an extra pole but its contribution may be discounted in the limit 
x, x' — > — oo, because it will be exponentially damped. Thus, the two integrals reduce to 

dk f dk' i i _i( fcx+fc v) ( 2ik \ / 2ik' 



2tt J 2tv u 2 _ k 2 _ 4 + ie u 2 _ k /2 _ 4 + ie \ik-2aj \ik' - 2a 

and hence to, 



9 B^'fe 1 ' (3 - 4) 



where = Vw 2 — 4. 



Combining the results contained in (|3.3| ) and ([Tj), the contribution (|3.1|) becomes 



ltan 2 a7r J 2n \ cosh 9 — sin a7r cosh ^ + sin a7r 

(3.5) 

where the integral over cj is to be understood as a positive energy integral, and it is 
convenient to set k = 2sinh# in the last terms of the integrand. The overall factor of a 
should be interpreted modulo unity. Summarising, the correction to the reflection factor 
from the boundary piece is 

^^ sinh e(—l _i ). (3.6) 

2 tan aix \ cosh 6 — sin an cosh + sin an J 



The evaluation of (|3.2| ) proves to be more testing. The integral over t" is straightforward 
again and provides a delta function allowing one of the non-loop energy integrals (u or u/) 
to be performed immediately. The loop integral is 



I dudk i (l + K.Jk")e~ 

J 2tt 2tt u" 2 -k" 2 -4 + ie V 



2ik x 



which is logarithmically divergent. However, a minimal subtraction of the term without the 
exponential dependence on a position coordinate leaves a finite integral. This subtraction 
is natural, coinciding with what would be achieved by normal-ordering in the usual bulk 
theory on the full line, as has been emphasised by Kim [jTB| . What remains is: 



. „ f du> dk 
2^ 

duj" dk" i 



J_ 00 J 2% 2% iv 2 -k 2 -4 + ie V K 1 J 

- Mfc") z~ 2tk " x " ( 3 - 7 ) 

ik'(x"-x') _|_ jQ a (J i '^ e -' ik ' i x '+ x ") 



2n 2n oj" 2 - A;" 2 - 4 + it 
du>' dk' i 



2tv 2tv u' 2 - k' 2 - 4 + ie 

The integration over x" is achieved using the device 

o 

dx" e {lX+p)x " = - 



A — ip 



where p is a positive constant which will be taken to zero at the end of the calculation. 
This introduces a collection of poles for the k" integration which need to be dealt with 
separately. 

The pole contributions after integrating out uj" are 

% fdudkdk' e _ Mt _ t > )e _ i{kx _ k > x > y i i 



2j2n2n2n~ w 2 _ jfc 2 _ 4 + ie u' 2 - k' 2 - 4 + it 



dk " ' fC a (k") 



2n Vk" 2 + 4 



jfc - 2k" + k' -ip k- 2k" - k' - ip 



-k - 2k" + k' -ip -k- 2k" - k' - ip 

and the k" integrations may be done by collapsing the contour into the upper half-plane and 
onto the branch cut running from k" = 2i to infinity along the imaginary axis, avoiding 
all the poles (the 'zp' effectively indicates the sense in which the k" contours must be 
interpreted), with the possible exception of the pole in JC a (k"). If a > this pole is 
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avoided automatically; if a < 0, it cannot be avoided but its residue integrated over k and 
k' yields exponentially decreasing terms as x,x' — > — oo. 

The integrals along the cut are all reducible to integrals of the following form 



dy 



1 



1 



1 



7T 



'y 



2-4 2/ + 2C 



tan" 

'1-C 2 V 2 



1 + C 



(3. 



Thus, for example, with a little manipulation, 



dk" 



2n Vk" 2 + 4 



K a {k>') 



1 



1 



k - 2k" + k' -ip n 



dy 



l 



k + k' - 2iy 



dy- 



1 



2 



2tv J 2 a ^J y 2_ A y y + i(k + k>)/2 
which can be evaluated using ( |3.8|) , to yield: 



y + 2(7 



i 

2^ 



2JC a ( 



k+k' 



/ 4+(k + k') 2 /4 



TV 

2 " tan " 



l4 + i(k + k') 
4-i(k + k') 



vT 



7T 

2 



tan 



l + <7 



1 



(3.9) 



The k, k! integrals can be performed next and, up to vanishingly small factors as before, 
one obtains as far as (|3.9|) is concerned a relatively simple result, 



i 

2r~v~ 



+ P V2 



an 



2 sin aw 



(3.10) 



The other three pole pieces may be treated in the same manner yielding contributions 
similar to ( |3.9|) , except that k + k' is replaced by one of k — k', —k + k' and —k — k', 
in turn, leading to corresponding differences in the analogues of ( |3.9p . Assembling all 
four contributions, and including the numerical factors from the integrations leads to the 
following correction to the reflection factor: 



i(3 2 C 2 



JCcr(k) sinh6* 



1 



1 



1 



2 Vcosh0 + 1 



+ 



a 



cosh 9 
1 



1 



sin 2 an \ cosh 9 — sin an cosh 9 + sin an 



(3.11) 
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In ( 3.11 ), the overall factor a is restricted to the range 0<a<l/2, in order that a be 
positive. 

The combinatorial factors C\ and C<i are equal since they arise from four-point inter- 
actions which are in essence identical, differing only by constant factors. Moreover, the 
common factor is 1/2, as is easily checked using, for example, the functional integral form 
of the Green's function. 



Finally, notice that the two contributions from the boundary summarised in (|3.6|) and 
Q3.ll ) combine neatly (within the restricted a range) to give the result: 



SK a (k) 



JC a (k) sinh6> 



1 



1 



cosh 9 + 1 cosh 9 
1 



+ 2a 



1 



cosh 9 — sin air cosh 9 + sin an 



(3.12) 



For the special case a = 1/2, (|3.12| ) collapses to (|L9|), as expected. 



4. Comparison with Ghoshal's formula 

To make the comparison with Ghoshal's formula up to 0(/3 2 ) the expansion of ( |1.5| ) is 
required. It is: 

/C(0) ~/C o (0) ^1 - ^sinh0 , (4.1) 

where 



1 1 

G{9) 



cosh 9 + 1 cosh 9 

ei e 1 



+ 
+ 



cosh 9 + sin(eo7r/2) cosh 9 — sin(eo7r/2) 
fi fi 



(4.2) 



cosh 9 + sin(/o7f/2) cosh 9 — sin(/o7r/2) ' 

and E ~ eo + ei/3 2 /47r, F ~ fo + fi(3 2 /4n. Comparison with the classical reflection factor 
JCq(9), eq (|1.13|) reveals that eo = Oq + oi, fo = ao — oi; comparison with eq ( |3.12|) reveals 



further that when ao = a± = a, 

P 2 



E ~ a 2 

V 2tt 

However, since fo vanishes, nothing can be concluded to this order concerning F. Never- 
theless, it is satisfying that the correction (|3.12| ) is of the same general form as (£ 
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Given the coupling constant always appears to enter in the combination B(/3), it is 
tempting to suppose on the basis of the above that 

£(a,/3) = 2a(l-|Y (4.3) 



On the other hand, consider the weak-strong coupling transformation applied to (|1.5|) . 
The numerator of (|1.5| ) is invariant under a duality transformation and, if it is assumed 
there is another theory within the same class whose boundary condition implies the dual 
reflection factor, there must be an a* such that 

either E(a,4n/P) = ±E(a*,/3) mod 4 or E(a,4^/P) = ±F{a*,j3) mod 4 . (4.4) 

This is a plausible assumption given the invariance of the two-particle S-matrix under the 
duality transformation, and given each boundary condition within the special class ao = a\ 
preserves the symmetry <fi — > — 0. The first possibility within ( [Oj ) cannot work given the 
form of (|4.3| ) — there is simply no suitable a*. On the other hand, the second alternative 
is possible and suggests 

B 1 

F(a,(3) = (l-2a)— and a* = - - a. (4.5) 

With F given by ( |4.5|) , there would be no change to the perturbation theory result ( |3.12j ) 
up to the order given. On the other hand, the next order should provide a test of the 
duality argument. 

However, the correct formulae cannot be so simple. Recall that a Dirichlet condition 
0(0, t) = also belongs to the special class, on allowing a to be pure imaginary and infinite. 
As before, the perturbative result may be obtained from (|4.2|) . On taking the appropriate 
limit, all that remains is the first part, in perfect agreement with Ghoshal's formula for 
the Dirichlet case, which is self-dual and given by 

(2-fl/2)(l + fl/2) 

A~D = 



(1) 

Within the special cases, it would be obtained by taking E — > zoo and F — > mod 4. 
Clearly, the correct formula would not be obtained using (fO|) and ( |4.5|) : ( f4.3|) is satisfac- 
tory but (|4.5| ) is not. 
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Moreover, the special coupling constant for which B = —2, the 'free-fermion' point of the 
sine-Gordon modeljl has been analysed before |H| and relations have been written down 
between the parameters r\ and d occurring in Ghoshal's formula at this value. Explicitly, 



using the notation of |T9| but in the normalisations of this paper, the ratio of the boundary 



to bulk coupling 7, and the parameter 0o are given by 

. , 1 / COSao7T COSOi7T 

7 = 4 cos ao7r cos ai7r, cosh (f>o = - \ h 



2 \cosai7r cosao7r 
In terms of these, 

k = (1 - 7cosh</> +7 2 /4)" 1/2 , cos£ = k(l - 7cosh0 o /2), 

and the two equations 

cos^coszi? = — (cos£)/fc, cos 2 ?] + cos 2 z?? = 1 + 1/fc 2 , 

are solved by 

cos 77 = ± cos(ao + ai)iv, cos = ± cos(ai — ao)7r. (4-6) 
Hence, it is simple to deduce that at the free fermion point, and within the cases ao = ai, 

^- = E = 4a ^=F = mod4. 

7T 7T 

Again, the formula (|4.3|) is adequate but the formula (|4.5| ) is not. 

A more general analysis might proceed along the following lines. Instead of (|4.3|) and 
(P~5|), one could write (q = /3 2 /4tt) 

B(a >5 ) = 2^ F(„, S )=2!^, (4.7) 
1 + q 1 + q 

in which case, the duality transformation would require, 

e(a,l/q) = ±f(a*,q) mod 2 ( J f(a,l/q) = ±e(a*,q) mod 2(1 + g). 



Q 

Then, taking e(a, q) = a would be a consistent choice (as suggested by the known facts 
concerning e(a,g)), provided 

a = ei/(e 2 /(o, 1/q) mod 2(1 + g), g) mod 2 f^^J > ( 4 -§) 



It is also a special point of the sinh-Gordon model in the sense that the S-matrix (1.1) is unity 
there. 
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where e\ and e 2 are choices of sign. However, it has not yet proved possible to obtain a 
satisfactory conjecture based on this idea. For example, it is straightforward to see that 
choices of /(a, q) of the form 

A{q) + aB{q) 
na ' Q) C(q) + aD(q) 

are inadequate once all the constraints are taken into account. 

Alternatively, one might return to the first of ( [4.4j ) and seek a formula for E(a, q) of the 
form 

, , = n (a) + m(a)g + n 2 (a)g 2 j + n k (a)q k 

[a ' q) d (a) + d 1 (a)q + d 2 (a)q2 + --- + d k (a)q k ' 1 } 

where the coefficients of the polynomials are constrained by the data described above 
and by the duality requirement. Neither linear nor quadratic polynomials appear to be 
satisfactory, while cubics and polynomials of higher degree are insufficiently constrained. 



5. Discussion 

Clearly, there is much left undone. Ghoshal's formula has been verified to lowest order 
and it has proved possible to find a compact parametrisation for the free-fermion point. 
However, it has not yet allowed a consistent conjecture for the general parameter depen- 
dence even within the special class of models considered here. In the general case, ao 7^ 01, 
the calculations are more intricate for the reasons stated before. 

It is intriguing to wonder about the supersymmetric sinh-Gordon theory. It has been 
noted [[20] that only a few boundary conditions are compatible with both supersymmetry 
and integrability. Moreover, they belong to the special class considered here, at least as 
far as the bosonic part is concerned (and in the normalisation of this paper correspond to 
o"o = o\ = ±1). It would be interesting to understand the effect of the fermions on the 
duality transformation, since one would guess the two allowed points are either self-dual, 
or dual partners. 

If and how duality works for the other affine Toda theories with a boundary will be inter- 
esting to discover. For the ade series of cases the S-matrix is self-dual, but the parameters 
which may be introduced at the boundary are highly constrained [H|7j , and known merely 
to be choices of sign. Presumably, the duality transformation could relate the possibilities 
in pairs, although some may be self-dual, and the Neumann condition is likely, bearing in 
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mind the result for sinh- Gordon, to be related to the choice (+) for all boundary parame- 
ters. The other affine Toda models are arranged in pairs related by duality and it is known 
that for most of them there are one or possibly two extra parameters permitted at the 
boundary. It will certainly be interesting to understand how duality will work for these 
cases. Up to now, some of the classical backgrounds have been found |f)lj21| but there has 
been no attempt to calculate the reflection factors perturbatively, except for the Neumann 
condition. 

Finally, one could take a different point of view. There remains the (unwelcome?) 
possibility that the weak-strong coupling transformation does not relate pairs of boundary 
conditions at all, that 



is the correct parametrisation, that the Dirichlet condition and the condition with a = 
are self-dual, while all the rest within the given class have reflection factors whose large (3 
limit is -1/(1) 2 . 
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